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Abstract 



To explore possibilities of avoiding coincidence problem in f(R) 
gravity we consider models in Einstein conformal frame which are 
equivalent to Einstein gravity with a minimally coupled scalar field. 
As the conformal factor determines the coupling term and hence the 
interaction between matter and dark energy, the function f(R) can 
in principle be determined by choosing an appropriate function for 
the deceleration parameter only. Possible behavior of f(R) to avoid 
coincidence problem are investigated in two such cases. 
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1 Introduction 

It is now more than a decade when the first observational evidence indicated 
that our Universe might have been expanding with acceleration [1J. With a 
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growing number of data favoring such evidences over past few years such 
doubt has now turned into certainty that the Universe has entered an accel- 
erated phase of expansion from a decelerated phase very recently[2j[3J. This 
observational fact however needs a solid theoretical support which appears to 
be a major challenge in cosmology today Efforts to resolve this matter have 
been made so far in two possible directions. One approach is to introduce 
a new term in the energy-momentum side of the Einstein's equation which 
will act as a source of gravity with repulsive effect and initiate acceleration 
in the course of expansion. Such a source has been termed in literature as 
'Dark Energy' due to its unknown nature. It is usually parameterized by an 
equation of state of the form u^e — p/p- From Friedmann equation one can 
see that a value of ude < — 1/3 is required for accelerated cosmic expansion. 
The most natural and simplest choice as a candidate for dark energy could 
have been the Cosmological Constant A with equation of state ude = — 1- 
But if one believes that vacuum energy is the origin of it then one fails to find 
any mechanism to obtain a value of A that is 120 orders of magnitude less 
than the theoretical prediction to be consistent with observation. Attempts 
were made by introducing the concept of varying A to explain observation but 
with same equation of state ujde = — 1 [4J . Possibilities of evolving code were 
also explored in many dynamical dark energy models. Primary candidates in 
this category are scalar field models such as Quintessence [5] and K-essence[6]. 
A common example to quintessence is the energy of a slowly evolving scalar 
field that has not yet reached the minimum of its potential V((f>), similar to 
the inflaton field used to describe the inflationary phase of the Universe. In 
quintessence models the parameter range is —1 < ujde < —1/3, and the dark 
energy density decreases with a scale factor a(t) as pde oc a~ 3 ( 1+UJDE ^[7\. A 
specific exotic form of dark energy denoted phantom energy, with ujde < — 1, 
has also been proposed[5]. It possesses peculiar properties, such as the vi- 
olation of the energy conditions and an infinitely increasing energy density. 
However, recent fits to supernovae, cosmic microwave background radiation 
(CMBR) and weak gravitational lensing data indicate that an evolving equa- 
tion of state crossing the phantom divide, is mildly favored, and several mod- 
els have been proposed in the literature [9J. In particular, models considering 
a redshift dependent equation of state, possibly provide better fits to the 
most recent and reliable SN la supernovae Gold dataset. A few compre- 
hensive and recent review articles on dark energy may be of interest in this 
regard [TO]. 

The alternative way to explain cosmic acceleration is to anticipate that 
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large scale dynamics of the Universe is not governed by Einstein's equa- 
tions. In a classical generalization of general relativity, the Ricci scalar R in 
Einstein-Hilbert action is replaced by a more general function f(R) which 
modifies the geometry-side of the Einstein's equation[llJ. Earlier interest in 
f(R) theories was motivated by inflationary scenarios as for instance, in the 
Starobinsky model[T2]. Recently, a number of models in f(R) modified theo- 
ries of gravity has been verified in an attempt to explain the late-time acceler- 
ated expansion of the Universe. In particular, it has been shown that cosmic 
acceleration can indeed be explained in the context of f(R) gravity [13], and 
the conditions of viable cosmological models have also been derived[14j. In 
the context of local tests in the Solar System regime, severe weak field con- 
straints however, seem to rule out most of the models proposed so far [To], 
although viable models do exist[16]. Some recent reviews on f(R) gravity 
may be of interest in this regard [T7]. 

There exist many cosmological models in 'modified gravity' as well as 
'dark energy models' that can successfully explain the observed late-time 
acceleration. However, there is no strong reason to choose one as better 
than the other. Not only that, there remain many unsolved issues related to 
all these approaches. One important issue among those is the Coincidence 
Problem[18j. If we accept that the observed acceleration is being caused by 
some 'dark energy' it is striking to find out that in our present Universe the 
density of dark energy and dark matter are of the same order. In fact any 
cosmological model successful in explaining late-time acceleration will neces- 
sarily suggest that dark energy density has started dominating over matter 
energy very recently. Now the question is why is it happening now? Is it a 
mere coincidence or is there some deep underlying reason behind it? This 
Coincidence Problem has been addressed in literature extensively [12]. One 
way to get rid of this is to design a model in which the ratio of matter den- 
sity to dark energy density settles down with evolution to its present day 
value quickly and remains stationary for a long period until today so that 
the coincidence problem does not arise. However such a model needs inter- 
action between dark energy and matter. Attempts to realize the coincidence 
problem as a consequence of interaction between matter sector and dark en- 
ergy are not new[20]. The suggested interaction terms in all such attempts 
are however phenomenological and are not derived from any fundamental 
theory. In this present note we address the coincidence problem in the con- 
text of f(R) gravity models in Einstein conformal frame. The conformally 
rescaled metric in such a model has a scalar field partner that act as another 
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dynamical variable of the vacuum sector. Understandably this scalar field 
is not a kind of matter field and is actually given in terms of f{R)- There 
exists a coupling between this scalar degree of freedom and the matter sector 
induced by this conformal transformation and the interaction between mat- 
ter and dark energy opens up the possibility of addressing the coincidence 
problem. The field equations look very much like those of quintessence mod- 
els but with a difference. The interaction term is given by the conformal 
factor and hence the coupling constant is universal in a sense that the scalar 
field couples with same strength to all kinds of matter fields. In a recent 
work Bisabr[21j claims to have derived necessary conditions for alleviating 
coincidence problem in some f(R) gravity models in Einstein frame. In these 
models the functional forms of f(R) are chosen by hand and the parameters 
in it are fixed by putting constraints to make the models free from coinci- 
dence problem and produce the present day observed value of deceleration 
parameter. We noticed however that the choice of f(R) fixes the dynamics 
of the model and although it produces desired value of present day decelera- 
tion parameter correctly; its time evolution appears to be not in agreement 
with observation. We suggest in this paper a different approach to overcome 
this difficulty We take up two models by choosing the functional forms of 
deceleration parameter from those available in literature and look for a suit- 
able f(R) that can alleviate the coincidence problem. In the first model the 
possible asymptotic behavior of f(R) that can solve the problem is shown. 
The second model however clearly rules out the existence of any well behaved 
f(R) for this purpose. 

2 Action and field equations 

In the f(R) modified theories of gravity, the standard Einstein-Hilbert action 
is replaced by an arbitrary function of Ricci scalar R and is given by 



where C m is the matter Lagrangian density and ip denotes all matter fields 
collectively. Using metric formalism, by varying the above action with respect 




(1) 



where k = 8nG. The matter action SM{g tiV ,i>) is defined as 
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to g^ v ', one can arrive at the field equations 



1 



where /' = df/dR. The matter stress-energy tensor is defined as 



rp(m) 



fJ,V 



2 S(^C r , 



(2) 



(3) 



The f(R) gravity may be written as a scalar-tensor theory, by introducing a 
Legendre transformation {R, /} — > {</>, U} defined as 



( j ) = f(R),u^) = R^)f , -f(RW)- 



(4) 



In this representation the field equations of f(R) gravity can be derived from 
a Brans-Dicke type action given by 



(5) 



This is the Jordan frame representation of the action of a Brans-Dicke theory 
with the parameter to — 0. As usual in Brans-Dicke theory and more gen- 
eral scalar-tensor theories, one can perform a canonical transformation and 
rewrite the action 03) in what is called Einstein frame. Rescaling the metric 

as 

-> g^v = f'(R)g^u (6) 

and redefining 0—7-0 with 




2k 



(7) 



the scalar-tensor theory can be mapped into the Einstein frame in which 
the 'new' scalar field couples minimally to the Ricci curvature and has 
canonical kinetic energy as described by the gravitational action 



S 



d 4 xy/—g 



|: - l -d^dJ - v{$) 



+ S M (e- 2 ^r u ^)- (8) 



The self- interacting potential V(<f>) is given by 

Rf'(R) - f(R) 



V(0) 



2kf'i(R) 



(9) 
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Apparently, a coupling of the scalar field (f) with the matter sector sector 

is now induced. The strength of this coupling (3 = y| is however, fixed 

and is same for all matter fields. Taking g^ and <p as two independent field 
variables, the variations of the action Qyield the following field equations 
respectively 

g, v = Hf^ + fp (10) 

= -pVkf- (ii) 



where 



T% = V M 0V,0-^V>V a 0-^(0)^ (12) 



~ -2 5S M (g^^) 



are stress-tensors for the scalar field and the matter field. It is important 
to note that T^ u and T™ v are not separately conserved and hence satisfy the 
following equation 

y^frn = _v^T / t = (3VkV u f m . (14) 

The trace of matter stress-tensor T m is obtained by taking trace of the field 
equation ffTOl 

f m = V a 4>Vj + 4V($) - R/k. (15) 

Hereafter in the rest of this article we shall stop using tilde overhead for 
convenient reading. 

For a spatially flat homogeneous isotropic Universe described by the 
Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (dr 2 + r 2 d6 2 + r 2 sin 2 6d(f) 2 ) (16) 

assuming matter in the form of pressureless dust the field equations f fTOj) 
reduce to 

3H 2 = k( Pm + P(f) ) (17) 
and 2H + 3H 2 = -ku^. (18) 
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Whereas field equation (fill) reduces to 

4> + 3H<j)+^- = -pVkp m . (19) 

d(p 

Here an over dot denotes derivative with respect to time t. The character p 
denotes the energy density of matter or other fluids. Throughout this article 
suffix m to a quantity denotes its matter part while the suffix denotes the 
'dark energy' part of it. H is the usual Hubble parameter defined as H = a/a 
and Utf) = is the equation of state parameter of the scalar field 0. The 
density p^ and pressure p^ of the scalar field are defined as 

= \tf + v{4>) 

and ptf, = ^0 2 -^O) 

Due to the coupling between the scalar field and matter field, the two 
components interact and hence are not separately conserved. From (ITT)) . 
(ITS)) and (TT9")) the conservation equations for matter and the scalar field follow 
easily to give 

p m + 3Hp m = Q (20) 
and p^ + 3H(l + u^)p ( p = -Q (21) 

where the interaction term 

Q = f3Vk<p Pm (22) 

This implies that for > 0, Q and hence p m is positive and energy flows 
from dark energy to dark matter and for < energy is transferred in the 
opposite direction. It is interesting to note that the interaction term (122)) 
looks similar to some of the phenomenological coupling terms suggested in 
literature [20]. It vanishes when = constant i.e. when f(R) is linear. A van- 
ishing Q implies that matter and dark energy remain separately conserved. 



3 Outline Of The Present Work 

To avoid the coincidence problem, matter and dark energy must scale each 
other over a considerably long period of time during the later stage of evolu- 
tion of the Universe. In other words, the ratio of two energy densities r = — 

P4> 
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remains constant in spite of their different rates of time evolution. Thus an 
interaction between the matter and dark energy becomes necessary in these 
kind of models. An appropriate choice of the interaction term may lead to- 
wards this goal but usually such a choice is phenomenological and does not 
come from any fundamental theory. However, in f(R) gravity models in the 
Einstein frame, the form of Q is determined by the transformation itself and 
one does not have the liberty of tailoring the interaction. Thus to achieve a 
stationary value for r and hence to avoid the coincidence problem, we can 
look for an appropriate form of f(R) instead. 

The set of three independent equations (ITT)) . (IT8l) and (TT9T) contains four 
unknowns viz: H, p m , and V(<p). Choice of any one of these parameters (or 
a derived quantity like the deceleration parameter q) shall in principle suffice 
to determine the functional form of all other quantities including f(R) and 
Q and one may check how well the model alleviates the coincidence problem 
fitting adequately with observations. In contrast, one can also start with a 
choice of f(R) and solve for the other parameters of the model using the set 
of equations presented in the previous section and check with observations. 
However, due to complexities involved in the field equations this approach is 
difficult to work out. In the present work we take the former approach i.e 
search for a function f(R) for a given form of q such that the ratio of energy 
densities r settles at a stationary value in course of evolution. 

The Hubble parameter H{t) = a/a and the deceleration parameter q(t) = 
(affU are related as functions of the redshift z(= —1 + by the equation, 



H(z) = H exp 



f Z {l + q{z')]dln{l + z') 
Jo 



(23) 



where the subscript refers to the present day values of the variables. Also 
note that by definition, 

H = -(l + q)H 2 . (24) 
Finally, the Ricci scalar R can be expressed in terms of q and H as 

R = Q(l-q)H 2 . (25) 
On the other hand, let us consider the time evolution of the ratio r = — , 

Pm p(j> /nc\ 

r = r— . (26) 

P4> P4> 



Using equations Q20J), (J2TJ) and (T22J) in (|26j) we obtain 
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r = "iHrujtf, + f3\fk(j)r(l + r) 
Application of ff24|) in equation (fl8|) leads to 



(27) 



2g- 1 



(28) 



Substituting (1281) in (1271) and using equation (ITTj) we arrive at 



r 



(2g- l)H + (3Vk<j). 



(29) 



r(l + r) 



We can show that the observed equality in the order of magnitude of the 
dark energy density and matter density of the Universe in recent epochs is 
not merely a coincidence if we demand that the ratio of the densities r either 
remain constant or vary extremely slowly as compared to the scale factor a 
over a long period of time. This means to analyze recent behavior of the 
Universe one can safely assume f = and therefore from (129]) we can write 



where a prime denotes derivative with respect to R. 

Thus for a given q(z) one can easily determine the evolution of the hubble 
parameter H and the Ricci scalar R as functions of the redshift z, using equa- 
tions (1231) and (125]) . These functions can then be used to solve equation (13TT) 
and a functional form of f{R), valid in the r = regime, can in principle be 
determined . In the subsequent sections we make attempts to determine f(R) 
for two different parametric choices of q(z) following the above prescription 
to alleviate coincidence problem. 



(2q - l)H = -f3Vk<p. 
Replacing in (1301) by using equation flTJ) gives 



(30) 




2(2q-l)H 



(31) 



4 Model I 



The simplest parametrization of the deceleration parameter q in terms of the 
redshift z available in literature is the linear two parameter relation given 
by[22][23] 
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q(z) =q + q'z (32) 

where q = q\ z= o is the present value of the deceleration parameter and 
q' = \ z =o- For a spatially flat Robertson- Walker Universe the best fit to 
the pair of parameters is (go, q') = (-0.73, 1.5) and the transition redshift for 
these set of values is zt = 0.49 [23]. Note that for z > 1.154, the deceleration 
parameter q > 1. Such a parametrization is therefore justified only at low 
redshifts over a small range < z < 1.15. Assuming q(z) as in ( 1321) . equations 
023]) and ([25]) give 



H(z) = H exp[g^](l + z) {1+q °- q,) (33) 
and R{z) = 6H%exp[2q'z}{l - q - q'z){l + z )^ 1+9> -^ (34) 

Integrating (I3T]) and by using (!32l and (133]) . we get 

| = ^exp[4g'z](l + z)^'- 2 )g (35 ) 

where f' Q = ^| 2=0 . Note that since the functional form of q(z) in this model 
is valid over a small range < z < 1.15, therefore the functional behavior 
of f{R) as depicted in the above equation is only asymptotic and may not 
indicate its true behavior for z » 1. 

The complexity involved in the expression on the r.h.s of (135]) does not 
permit an analytical solution for f(z). However, since the expressions are 
true only for small z we can always look for an approximation. Note that 
the expressions of H(z), R(z) and j- given by ( 1331 . f l34"]) and fl35]) contain 
an exponential term each. Retaining the terms up to the sixth power of 
z in the series expansion of the exponentials and using them in f|3B]) an 
analytical form of f(z) can be obtained. The analytical expression, however, 
is not being mentioned here for economy of space. Instead variation of f(z) 
versus z under the above approximation is plotted in figJU That such an 
approximation holds good in our range of interest (0 < z < 1.15) is also 
clear from figffia)and (b) where we find that the curves for H(z) and R(z) 
drawn with and without the approximation appears almost indistinguishable 
at low redshifts. A parametric plot of R and / with z as parameter, is then 
plotted in fig]]Jd)to obtain the nature of variation of f{R) against R. 

Thus in order to address the coincidence problem in this model the func- 
tion f(R) should be such that, at low redshift it behaves like the curve given 
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in figJT^d). The plot shows that the function f(R) which is initially linear in 
nature for z > 0.5 implying no interaction with matter sector, attains a non 
linear character in the very recent past z < 0.5 facilitating the interaction 
and hence a transition from deceleration to an accelerated phase of expansion 
of the Universe. It is always possible to suggest an analytic function that can 
match the behavior of f(R) as plotted in fig. 1(d) but we must remember that 
any such choice that fits the curve does not guarantee r — 0. For example 
the function 

f(R) = XiR + X 2 R n (36) 

when plotted against R appears to behave similarly (shown in figE]). 

The behavior of r for the above choice of f(R) however deviates from 
desirable features. For this reason we refrain from suggesting any f(R) but 
stress on the fact that such a function can always be found in principle that 
not only fits the curve but also guarantees r = 0, thereby alleviating the 
coincidence problem. 



5 Model II 



Another parametrization of q(z) common in literature is of the form|24j 



q(z) 



1 glf + g2 

2 (1 + zY 



(37) 



where q± and q% are constants with (| 
deceleration parameter and \ 



.., , q 2 ) giving the present value of the 
| 2 =o = Qi — 2g 2 - For a spatially flat Robertson- 
Walker Universe the best fit to the pair of parameters is (qi, q 2 ) = (1.47, 
-1.46) and the transition redshift for these set of values is zt = 0.36[25]. For 
high redshifts z ^> 1, the deceleration parameter q remains fixed at | in this 
model. 

Proceeding in exactly the same manner as in section (j4]) we get 



H(z) 



R{z) = 6# 2 (1 + z) 



and 



1 
2 

dz 



H Q (1 + z)* exp 
qiz + g 2 



02 q\z^ - q 2 
2 2(1 + z) 2 



(1 + 2 

/oexp 



exp 



12 + 



(i + zY 



-2{2q lZ + qi + q 2 ) 
(l + zf 



dR 

dz 



(38) 
(39) 
(40) 
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(a) (b) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 5 6 8 9 10 

2 R 



Figure 1: The plots (a) and (b) show the variation of actual and approximate 
expressions of H(z) and R(z) vs. z in red and blue respectively for model I 
with q = —0.73 and q' = 1.5. The curves appear to match exactly at low 
redshifts. Variation of f(z) vs. z and f(R) vs. R are shown in (c) and (d) 
respectively with f' Q = 2 and H Q — 1 units. 
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(a) (b) 




0.5 1 1.5 2 



Figure 3: The plots (a) and (b) show the variation of actual and approximate 
expressions of H(z) and R{z) vs. z in red and blue respectively for model II 
with qi = 1.47 and q<i = —1.46. The curves appear to match exactly over a 
long range of redshift. Variation of f(z) vs. z is shown in (c) for f = 2 and 
Ho = 1 units. 



Obtaining an analytical solution for f(z) from equation ( 1401) is difficult in 
this well. However, following similar approximation in the exponential 

terms of the above equations, as in the previous model, a solution for f(z) can 
be easily obtained and the nature of the function is shown in figE] along with 
the variations of H (z) and R(z) against z. The plots of H(z) and R(z) with 
and without the approximations are found to be almost indistinguishable over 
a long range of z justifying the approximation. Nevertheless, the parametric 
plot in figlU shows that the function f(R) in this case is double valued in the 
recent phase of evolution of the Universe. 

Therefore, for the particular parametrization of q(z) chosen in fl37|) . one 
cannot obtain a suitable well behaved function f{R), capable of alleviating 
the coincidence problem. This result in contrast with that obtained in[2T] 
can be understood if we notice that a prior choice of f{R) fixes the functional 
behavior of q and hence that of q' . In[2T] f(R) is chosen by hand and the 
present value of q is used to fix the parameters in the given f(R), but it fails to 
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Figure 4: Plot of f(R) vs. R for model II with q\ = 1.47, q2 = —1.46 and 
/o = 2 in H — 1 units. 



achieve the correct value of q'\ z =o in the process. Our approach on the other 
hand, is just opposite. We fix the function q that fits with observational data 
and look for a suitable f(R) that can alleviate the coincidence problem. We 
find that no such f(R) exists as long we restrict ourselves to q that behaves 
as given in (1371) . 

6 Conclusions 

A suitable kind of interaction between dark energy and dark matter can 
make the ratio of their densities possible to attain a stationary value during 
the course of evolution and consequently can help alleviating the coincidence 
problem. In Einstein gravity such interactions are phenomenological and 
needs fine tuning of the interaction term. On the other hand, in Einstein 
frame representation of f(R) gravity models, the interaction between the 
dark energy (scalar partner of the metric tensor) and dark matter is fixed by 
the conformal transformation. In this present article we address coincidence 
problem in the context of f(R) gravity in Einstein frame. We point out 
that fixing the functional form of any one of the four quantities viz. the 
Hubble parameter, the matter density, the scalar field or the potential (or 
any derived quantity like the deceleration parameter ) shall in principle suffice 
to determine the functional form of all other quantities including f(R). We 
choose two functional forms for deceleration parameter readily available in 
literature those fit observational data. The linear function q(z) in Model 
I determines the the functional behavior of f(R) for low redshift which is 
plotted in figJTJd). The plot shows that the function f(R) which is initially 
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linear in nature for z > 0.5 implying no interaction with matter sector, 
attains a non linear character in the very recent past z < 0.5 facilitating 
the interaction and hence a transition from deceleration to an accelerated 
phase of expansion of the Universe. In another attempt a popular non- 
linear function of q(z) is chosen for investigation. However the function f(R) 
obtained in this case is double valued at low redshifts and hence is not suitable 
from physical point of view. Thus our investigation rules out the possibility 
of such form of deceleration parameter in f(R) gravity if one is looking for a 
model free from coincidence problem. 
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